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ABSTRACT 



Context. The YohkohfHXT observations of Long Duration Events (LDEs) showed that the HXR emission (14-23 keV) is present for 
tens of minutes after flare maximum. Hence, some heating process is expected to exist during that time. The better energy resolution 
of RHESSI compared to HXT allow us to analyse LDEs in more comprehensive way. 

Aims. We selected three LDEs observed by RHESSI to answer the questions how long HXR emission can be present, where it is 
emitted, what is its nature and how much energy should be released to sustain the emission. 

Methods. We used RHESSI data to reconstruct images of the selected flares with an energy resolution as good as 1 keV. Next we 
estimated physical parameters of HXR sources through imaging spectroscopy. The physical parameters obtained were then used to 
calculate the energy balance of the observed sources. 

Results. We found that HXR thermal emission can be present for many hours after LDE flare maximum. The emission comes from 
large and hot loop-top sources. The total energy that must be released to sustain the emission of the sources is as high as 10^^' erg. 

Key words. Sun: flares - Sun: X-rays, gamma rays 



1. Introduction 

A long duration event (LDE) is a solar flare characterized by a 
slow decrease in soft X-ray (SXR) thermal emission. This de- 
crease may last from several hours to more than a day. Much 
insight into the nature of LDEs was made by ultraviolet and 
X-ray observations, during the Skylab, S MM (Solar Max imum 
Mission) and Yohkoh space missions (e.g. Sheelev et al., '1975; 
jKahler, 1977; Feldman et al., 1995; Harra-Mumion et al., 1998; 
ICzavkowska et all 119991: iShibasakil 120021: ilsobe et all 120021) . 
One of the most important conclusions is that without the con- 
tinuous energy input during the whole decay phase LDEs would 
decay much faster than it is observed. 

Loop-top sources (LTSs) are remarkable SXR and HXR 
(hard X-rays) features of solar flares seen close to a flare loop 
apex. They form before fla re maximum and in LDEs may las t 
many hours (e.g. Feldman et al.L 119951; iKolomariskil l2007ah . 
LTS s should b e located close to the primary energy release site 
(e.g.l Kopp & Pneuman, 1976; Shibata, 1999; Hirose et al., 2001; 
iKar lickv & Barta, 2006). Thus, they are a very promising tool 
for an analysis of an energy release during the decay phase. 
LTSs were first recorded on images taken from Skylab and are 
commonly present in the Yohkoh/SXT flare observations. Since 
the first observation it has become clear that the presence of an 
LTS during the whole flare decay-phase requires continuous en- 
ergy release and some restriction mechanism efficiently prevent- 
ing ou tflow of mass and energy from an LTS (see Vorpahl et al., 
Il977h . Without meeting these two requirements loop-top sources 
would rapidly lose energy by radiative and conductive processes 
and would vanish. This result was later con firmed by an anal- 
ysis based on RHESSI and Yohkoh data (e.g. Ijiang et all l2006t 
iKolomaiiski, 2007b). 

In the CSHKP (e.g ICarmichae l ITgel iHiravamaL [19741: 
iKopp & Pneumad. Il976l : Tsturroch Il966h model of a so- 



lar flare the energy release occurs above an LTS . This 
model and its numerous modifications ( Kruc ker et al.L l2008l 
and references within) successfully explain observational 
features of solar flares observed in different wavelengths 
during the impul sive phase. However , the authors inves- 



tigating LDEs (NakaiimaetalL 119981; lUchida et all 119991; 



Kundu & Grechnev, 2001; Philli ps et all l2005h report observa 
tions inconsistent with the CSHKP concept. The main problems 
they had were related to non-thermal velocities of 20-50 km/s, 
long times needed to keep the inflow of a magnetic field into a 
reconnection site and a lack of downflows or upflows. Recently 
Jiang et al. (2006) investigated heating and cooling processes of 
an LTS and found that there is a large amount of energy released 
during the decay phase of a short-duration solar flare. The au- 
thors suggest that in the case of an LDE the energy released dur- 
ing the decay phase may be even larger than the energy released 
during the impulsive phase. 

An analysis of energy release during the decay phase can 
give us precise constraints for numerical flare models. The most 
demanding for the models is long-lasting HXR emission of LTSs 
seen in LDE flares. Therefore such LTSs observed for many 
hours after a flare maximum are probably the most promising ob - 
servational feature to set up the constraints. The sources should 
be very close to the energy release site and, that is the most im- 
portant attribute, they put high requirements on the energy re- 
lease rate. If a long-lasting HXR source is thermal then it must be 
continuously heated, because the characteristic radiative cooling 
time of hot (above 10 MK) and dense (^ lO'" cm"-') plasma is 
X 1 hour If an HXR source is non-thermal then there should be 
a continuous acceleration of particles, to counteract fast thermal- 
ization of non-thermal electrons. Time-scale of the thermaliza- 
tion in the sense of momentum loss of electrons is about several 
seconds. 
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Hard X-ray LTSs were reported before with Yohkoh/HXT 
data, but with observations shorter than thermal soft X-ray 
LTSs. iMasuda et al. (1998) observed HXR sources up to 30 
minutes after flare m aximum, with typical sizes 1-2 arcmin. 
iHarra-Mumion et"aL I fl998 ) investigated two LDEs and con- 
cluded that HXR sources grow with time, with diameters x 
20 - 45 arcsec. In one of these flares the HXR source was visible 
up to 3 hours after flare maximum. Kolomanski (2007 a.b) anal- 
ysed three LDEs using Yohkoh data, finding that HXR sources 
during the decay phase were well correlated with SXR LTSs and 
were observed (in 14 - 23 keVrange) for about 50 minutes after 
flare maximum. ' Khan et akl (1200 6) reported the observation of 
thermal HXR source lasting more than 1 hour after flare maxi- 
mum. In all these cases the HXR emission was observed close to 
and above the SXR emission source. Although YohkohlVDiT data 
have provided useful information about LDEs, generally only L 
channel (14-23 keV) showed significant emission for any length 
of time. It was thus not possible to distinguish between the ther- 
mal andjioi>thermaljiato of an LTS. 

iGallagher et al] (l200l observed evolution of an LTS using 
RHESSI and TRACE data. The HXR emission was observed in 
the 12-25 keV range some 4 hours after the flare maximum and 
almost 1 1 hours in the 6-12 keV range. The LTS was large and 
its altitude increased at a speed that gradually declined from 10 
to 1.7 km s"'. The higher energy emission (12 - 25 keV) was 
located above the lower energy ones (6-12 keV and 3-6 keV) 
and all HXR emission was located above the tops of the loops 
obsei-ved in the EUV range {TRACE 195 A). 

So far, most investigations of LDEs made with RHESSI have 
been rest ricted to a few wide (s everal keV) energy intervals. 
Recently ISaint-Hilaire et al] (1200 9) presented observations of a 
coronal source recorded by RHESSI. The source was observed 
for about 12 hours in energies up to 10 keV. Emission of the the 
source was purely thermal with maximal temperature of about 
11 MK. Authors estimated a total energy (about 4 x 10-" erg) 
which had to be deposited in the source to sustain its long last- 
ing emission. 

Here we present the investigation of LDEs observed by 
RHESSI. The analysis is made using RHESSI images recon- 
structed in very narrow (1 keV) energy intervals. RHESSI is very 
useful for analysing weak sources even with 1 keV energy reso- 
lution due to its high spectral resolution. Thus, it gives an oppor- 
tunity to determine what is the nature of LTSs emission - ther- 
mal or non-thermal. Using the images we estimate LTSs physical 
parameters through imaging spectroscopy. Then the parameters 
are used to calculate the energy balance of the observed sources 
and to find out how effective are the energy release and heating 
processes at the decay phase of an LDE. 

2. Observations 

Three LDEs, well observed by RHESSI dLin et al.Ll2002l) . were 
selected for an analysis. We chose flares of significantly differ- 
ent power and with decay phases lasting more than 10 hours 
in GOES 1-8 A range. The flares occurred on: 30 July 2005 
(hereafter the LDEl), 22 August 2005 (LDE2), 25 January 2007 
(LDE3). Basic information about the LDEs is given in Table [1] 
Our analysis is based on RHESSI data (Reuven Ra maty High 
Energy Solar Spectroscopic Imager, iLin et al.l l2002i) supported 
with SoHO/EYT (Extreme UV Imaging Telescope on-board 
Solar and Heliospheric Observatory, iDelaboudiniere et"an . 
Il995h and GOES/SEM (Space Environment Monitor on- 
board Geostationary Operational Environmental Satellites, 
iDonnellv et al.L [T977h observations. 



2.1. On the RHESSI data 

RHESSI is a rotating Fou rier imager with n ine detectors made of 
pure germanium crystals (iLin et al I I2OOI . The detectors record 
energy and arrival time for each detected HXR photon. For 
strong solar flares the number of coming photons easily reaches 
10^ counts per second. As the number of cou nts rises the hfetime 
of a detector decreases. The pulse pile-up (ISmith et al.L l2002b 
is the main problem connected with this effect. It occurs when 
two or more photons arrive almost at the same time. In such a 
case the electronics of the detector recognize them as one pho- 
ton with energy equal to the sum of these photons. When number 
of incoming photons is very high the pulse pile-up may severely 
influence measured fluxes. 

The number of HXR photons decreases strongly with energy 
so problems with pulse pile-up occur mainly in the low-energy 
end of the range observed by RHESSI. Thus, to achieve the great 
dynamic range of the instrument, several systems to decrease 
the measured low-energy HXR photons were included. Among 
these systems there are attenuators which are the most impor- 
tant source of uncertainty. Influence of the attenuators on mea- 
sured fluxes is quite well understood at present, but in practice 
the attenuator cou ld seriously complicate the spectral analysis 
dSmith et al.Ll2002l) especially in the lowest energies. 

For each orbit, RHESSI passes through the radiation belts 
and South Atlantic Anomaly (SAA) which may have also a sig- 
nificant influence on the measured HXR fluxes. In our analysis 
we focused on the late phase of an LDE when the emission is ex- 
tremely weak, so some care is needed in the analysis. However, 
the image spectroscopy we performed has an advantage in com- 
parison to "standard" spectroscopy based on fluxes measured for 
the whole Sun. The background photons do not influence mod- 
ulation profile for the rotating Fourier imager. Thus, an anal- 
ysis of very weak fluxes can be done much easier than in the 
case of "standard" spectroscopy, especially for moments of pas- 
sage through radiation belts when the background changes very 
quickly with time. 

RHESSI light curves of the LDEs are shown in Figs. |3]|H] and 
fT2l The effect of radiation belts and SAA passages can be seen 
on the light curves. Usually the SAA influence is observed as a 
rise of the signal observed several minutes before and after the 
passage. When satellite is in the SAA detectors are switched off. 
The influence of the radiation belts is visible as a slow oscillation 
of the signal. 

3. Analysis 

3.1. Imaging 

It is a great difficulty to reconstruct a RHESSI image when an 
emitting source is weak. Usually RHESSI image reconstruction 
is performed for detectors Nos. 3-6, 8 and 9. Depending on the 
reconstruction algorithm and weights chosen this set of detectors 
gives a spatial resolution about 7-9 arcsec (lAschwanden et all 
I2OO2I) . However, we were not able to reconstruct any reliable 
image in a case of very late phase of LDE decay with the use of 
this standard set of grids. The weak signal was not the problem 
since taking integration times up to 8 minutes we were able to 
collect enough counts, i.e. about several thousand. 

The problem may be solved if we remember that when the 
source size is comparable to the resolution of a particular grid 
then the detector records very weak or no modulation of the sig- 
nal ( iHurfordet al.Ll2002h . The resulting image obtained for this 
grid is then dominated by noise. Such grid added to the set used 
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Table 1. Basic information about the selected long-duration flares. 



Flare 


Date 


Location 
Active Region' Coordinates^ 


SXR Flux Maximum^^ 
Time [UT] Class 


Duration^^ 
[h] 


Time Intervals Analysed 
Beginning [UT] (Duration [min]) 


LDEl 


30 July 2005 


10792 N10E59 


06:36 X1.3 


11 


10:20(0.07), 11:30 (0.13), 
11:55 (0.13), 12:40 (0.27), 
13:34 (0.53), 14:15 (1.00), 
15:55 (4.00), 16:26 (4.33) 


LDE2 


22 August 2005 


10798 S10W52 


01:34 M2.7 


11 


01:35 (0.13), 01:55 (0.27), 
02:45 (2.00), 03:24 (1.00), 
04:22(1.00), 04:46 (4.00), 
05:56 (8.00), 06:36 (2.00), 
07:29 (2.00), 08:08 (2.00) 


LDE3 


25 January 2007 


10940 S07E90 


07:15 C6.3 


17 


07:03 (1.50), 07:55 (0.13), 



08:12 (0.33), 08:32 (0.50), 
09:30 (1.50), 09:50 (1.50), 
10:10 (1.50), 11:00 (2.50), 
11:22 (3.00), 11:47 (4.00), 
12:30 (4.00), 12:56 (4.00) 



NOAA active region number; heliographic coordinates; according to GOES 1-8 A database 



for reconstruction introduces only noise. In the worst case we 
may not get convergence of the reconstruction algorithm. 

We tried to solve this problem in the following way. First, 
we reconstructed images for single grids and a wide field of 
view (8 arcsec pixel, image size 256x256 pixels) using the back- 
projection algorithm (Hurfordet al., 2002). Next, we used these 
images for determining which grid provides us with reliable im- 
ages of the source. An example is presented in the Fig. [T] For 
some time intervals there is only noise in the narrow grids, al- 
though count rates are high. Taking this into account we selected, 
for a given time interval, only those grids that showed a definite 
source in a single detector image. 

In the next step we used chosen set of grids to rec onstiTict 
images with PIXON algorithm ( Puett er & Yahill 1 1 999l and ref- 
erences within). Since integration times are rather long (min- 
utes) we used stacked modulations Q. Images obtained with this 
method were more "stable" than without it. "Stable" means that 
images in neighbouring time and energy intervals were similar. 
Moreover there were fewer images where we did not get conver- 
gence of the reconstruction algorithm. Energy resolution chosen 
for reconstruction was as narrow as possible, i.e. 1 keV. It en- 
abled to obtain reliable fits to observed spectra. 



3.2. Geometrical properties ■ 
and size 



LTS centroid location, altitude 



Having reconstructed set of images we were interested in ob- 
taining some properties of observed sources. In each image we 
determined sources using 50% of brightest pixel intensity iso- 
Une. In each case we observed sources related to coronal parts 
of flaring structure only (loop-top sources). There was no ob- 
served emission localised on the chromosphere level (footpoint 
sources) which is rather typical for the decay phase of an LDE. 

We determined position of loop-top sources by finding their 
centroids. The centroids were estimated as centres of gravity 
of the emission within the 50% intensity isoline (relative to 
the brightest pixel in the RHESSI image). The uncertainty of 
the source location estimation is caused by two agents. First 



' http://sprg.ssl.berkeley.edu/~tohban/nuggets/ 
?page=article&article_id=39 



agent is an uncertainty connected with spatial resolution of the 
instrument although using centroids we significantly improved 
the accuracy of source location estimation with regard to nom- 
inal. In our analysis this error was always less than 1.5 arcsec. 
More uncertainty was produced by the behaviour of analysed 
sources. We observed often a systematical trend in source lo- 
cation. Namely, the higher was the energy of the source the 
higher in the corona it was observed. Such systematical effect 
should be considered when estimating an actual source position. 
Thus, for a particular time interval we measured a spread of val- 
ues obtained for different energies and the spread was assumed 
as a good approximation of uncertainty in source position. We 
usually obtained values about 5 arcsec. This systematical trend, 
higher the energy h igher the altitude, is a common observational 
phenomenon (e.g. Gallagher et al., 2002; Sui & Holman, 2003t 
Liu et al., 2008; Saint-Hilair e et a l., 2009). Usually centroids of 
a source in energies below 10 keV are up to several arcsec below 
centroids in energies about 20 keV. 

To estimate the altitude of an LTS we need to know its po- 
sition in the RHESSI image and find the point on the solar pho- 
tosphere above which the source is situated (reference point). 
The source position was defined, as mentioned, by the position 
of its centroid. The position of reference point was taken from: 
locations of foot point sources observed by RHESSI at impulsive 
phase (LDEl), locations of flar e ribbons recorded by S oHO/EYT 
(LDE2), method described bv iRov & Datlowd (Il975h (LDE3). 
The method allow to estimate heliographic coordinates for flares 
behind the solar limb. We plot a position (on the solar disk) of 
an active region in which an analysed behind-the-limb flare oc- 
curred as a function of time. The position is taken from positions 
of all on-disk flares which occurred in that active region. Then 
we extrapolate the position vs. time plot behind the limb to get a 
position of our behind-the-limb flare. 

Since a moment of time for which position of the reference 
point was determined is different from times intervals analysed, 
all the positions were coiTected to account for the solar rotation. 
Finally, the LTS altitude was calculated as the distance between 
its centroid and the reference point. The altitudes obtained were 
corrected for the projection effect. Errors in the location of ref- 
erence point and the location of source centroid were included 
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Fig. 1. The 6 - 7 keV back-projected RHESSI 
images reconstructed for selected grids (Nos. 
5, 6 and 8) and three time intervals (the 30 July 
2005 flare). When a source size is comparable 
to the resolution of a particular grid then image 
obtained for this grid is dominated by noise. 
Field of view of each image is 1024x 1024 arc- 
sec. The Sun limb is marked be a black line. 



in the altitude errors. The relative error of altitude is about 10 
percent. 

The actual size of an LTS was determined using PIXON im- 
ages reconstructed with the use of all detectors in which we ob- 
served the source at a given time. The size of the source was 
estimated as the FWHM size. We did not observe important and 
systematical change of source size with energy. Thus, we de- 
cided to use PIXON images reconstructed in the energy interval 
7-8 keV for the estimation of a source size. 



3.3. Imaging spectroscopy of an LTS 

The main motivation for reconstructing images in narrow energy 
bands is the spectroscopy of the spatially resolved HXR sources. 
For this purpose we used OSPEXQ package. HXR continuum 
was fitted with thermal and, if it was observed, non-thermal com- 
ponents . We usually observe d also line features at 6.7 keV and 
8 keV ("Phillips e t al.L l2006l) . Physical parameters of observed 
sources were obtained from fits. Namely, for each time interval 
we obtained temperature and emission measure (thermal emis- 
sion). If a non-thermal component was observed then we fitted 
it with double power-law function. We treated a break energy 
and power law index above this energy as free parameters. The 
power law index below the break energy was fixed at constant 
value equal to zero. 



3.4. Energy balance 

As mentioned in the Introduction, the presence of HXR emis- 
sion from an LTS during the decay phase is evidence for energy 
release at that time. To calculate heating rate of an LTS we con- 
sidered its energy balance during the decay phase. Three major 
cooling processes where included into this balance: expansion, 
radiation and conduction: 



d£\ ld&\ 



where: 



(1) 



^ http://hesperia.gsfc.nasa.gov/ssw/packages/spex/doc/ 
pex_explanation.htm 



- & - 3NkT is thermal energy density, 

- ^ is the decrease of & per second estimated from tem- 
perature (T) and density (N) values, 

- ^ is the decrease due to the adiabatic expansion of 
plasma in a source, 

- Ec is the energy loss due to thermal conduction, 

- Er is the radiative loss, and 

- Efj is the heating rate or thermal energy release. 

The values of Ec, Er and Eh are in erg cm"^ s"'. We calcu- 
lated: 

-(fL = 5^^(f)' 

- Ec = 3.9 X \Q--''T^-^ l{Lr) where r is the LTS radius and L is 
loop semi-le ngth (assuming that e nergy is conducted to the 
chromosphe. lJakimiec et al.lll997h . and 
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- Er = N^(T) where (r) is the radiative loss function taken 
from lDere et al ] (l2009h . 

We took the altitude of an LTS above the photosphere 
(h) as an approximation for L in the expression for Ec- Of 
course, h is smaller than L, but they do not differ too much. 
From Equation ([TJ we calculated upper and lower limits for 
Eh- The upper limit, {EH)max, is calculated directly from 
Equa tion ([B i.e. wit h Eg d escribed by Spitzer thermal conduc- 
tion dJakimiec et all Il997h . However, it has been shown that 
the actual conduc t ivity m ay be smaller than Spitzer conductiv- 
ity dLuciani et al.L Il983h . A lower (the lowest possible) limit, 
(£■//)„„„, can therefore be obtained assuming Ec = 0, though this 
may not be physically realistic. The actual value of LTS heating 
rate must be contained between the upper and lower limit of Eh- 



3-5- Diagnostic diagram 

We used diagnostic (temperature-density, T - N) diagrams as 
an independent tool for analysing energy release. The method 
of investigating flare e volution by tert iperature-densi ty dia grams 
was put forwa r d by J j ^miec et al.l (Il986i 1 19871 [l992 ) and 
ISvlwester et al.l (Il993h . It was shown that temporal evolution of 
flares does not depend much on their morphology. Thus, the di- 
agrams allow us to study temporal changes of the energy release 
in flar es qualitatively (Jakimiec & Fludra, 1992). Jakimiec et al. 
(Il992h in their hydrodynamic modelling considered a single, 
symmetrical loop with constant cross-section and constant heat- 
ing rate before a flare. The authors took the following variation 
of the heating rate with time t in their calculations: 



7.5 



Enit) 



A - const for < f < fi 
Ae^~^ forf>fi 



(2) 



and studied two evolutionary tracks with r = s and t - 300 s 
for a flaring loop of semi-length L - 2x 10'* km. 

Non-flaring, steady loops having the same semi-length L but 
different temperature T and density lie along a straight line 
(the line of steady state loops, S-S line). Dependence between 
T an d for the S-S coronal loops is described by "scaling law" 
(Hosneretal., 1978): 



6.9 X 10"'^(A^L) 



1/2 



(3) 



0.5. A flare 



The inclination of the S-S line is ^ = Sof^ 
evolution on the diagnostic diagram can be divided into few 
stages (see Fig.|2]l. 



1. AB - Point A, lying on the S-S line, represents pre-flare 
conditions in a loop. Heating rate is constant, the loop is in 
steady-state. At the beginning of the flare energy is released 
abruptly causing increase of temperature (path AB). 

2. BC - Temperature is constant, electron density increases 
gradually due to chromospheric evaporation. Rate of energy 
release starts to decreases at point C. 

3. CD - Decreasing heating rate (point C) causes decrease of 
the temperature. However, temperature is still high enough 
to cause more chromospheric evaporation and increase of 
electron density (path CD). When released energy becomes 
insufficient the density in the loop begins to decrease (near 
the Une S-S, point D). 

Evolution of a flare during the decay phase depends on how 
fast the heating rate Eh decreases. In calculations two boundary 
scenarios observed during the decay phase may be present. 
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Fig. 2. Density-temperature diagnostic diagram illustrating the 
dependence of the slope on the characteristic decay time of the 
heating rate (t = s a nd t - 300 s cases are shown). The figure 
comes from the paper lj^imiec et aH (Il992h . 



- If Eh decreases very fast, with an e-folding decay time 
shorter than the thermodynamic decay time t,/, oc L/ -yTo 
(L - semilength of a loop. To - initial temperature at the top 
of a flaring loop) then the flare evolves along the OFF branch 
{Eh switched-ofF evolution). The evolution is characterized 
by a fast decrease of T caused by the loop cooling due to 
thermal conduction losses and radiative losses not balanced 
by the heating rate, which is too low. This evolution is rep- 
resented in Fig. |2]by a branch for Eh e-folding decay time 
T = Os. 

- If Eh decreases slowly, with r longer than T,h, the flaring 
loop will go through a sequence of steady-state configura- 
tions. In a r - diagram such a flare evolves along the QSS 
branch (quasi-steady-state evolution) which is parallel to the 
S-S line. During the QSS evolution T drops slower than for 
the OFF evolution because the thermal conduction and ra- 
diative losses are balanced by the heating rate which is high 
enough. The QSS evolution is shown in Fig.|2]by a branch 
for the Eh e-folding decay time t - 300 s. 

These results were obtained for a loop with a constant length. 
However, observations show that the height of flaring loops usu- 
ally increases during a flare evolution. If density also decreases 
(which is typical behaviour of solar flares) then QSS path will 
have inclination smaller then 0.5 (see Kolomanski, 2007b). In 
typical LDEs this altitude-density factor can lower the inclina- 
tion of QSS path by about 0.1 - 0.2. 

The time behaviour of the energy release during the de- 
cay phase of the analysed flares was investigated using simpli- 
fied diagnostic diagram (log T vs. log EM) obtained from the 
GOESISiEM data. Diagnostic diagrams obtained are shown in 
Figs.|7][II]and[l5l 

4. Results 

4-1- Ttie 30 July 2005 fiare 

RHESSI light curves of the LDEl are shown in Fig.|3] The flare 
emission could be analysed for eight intervals, given in Table [1] 
and|2] In the lower standard energy range 6-12 keV, RHESSI 
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Table 2. Basic parameters of the loop-top source of the 30 July 
2005 flare. 



06:00 08:00 10:00 12:00 14:00 
Start Time (50-Jul-05 05:00:02) 



16:00 



Fig. 3. RHESSI light curves for the 30 July 2005 flare. Vertical 
lines mark the boundaries of the satellite night (N) and the South 
Atlantic Anomaly (S) periods (solid and dashed lines mark the 
beginning and the end of the periods, respectively). The light 
curves were shifted vertically for for clarity. The intensity is 
given in relative units. 



recorded the emission for 10 hours during the decay phase un- 
til the beginning of the next flare. In the higher energy range 
12-25 keV the emission was detectable for 6 hours after flare 
maximum. 

Appearance of the flare and its evolution are shown in Fig. 2] 
The figure presents BIT 195 A images with a source observed by 
/?//£'5'5'/ overlaid on. As can be seen in EUV the LDEl consisted 
of several loops forming an arcade. In X-rays there was one loop- 
top source located just above the arcade of loops visible. 

Basic geometrical and physical properties of the loop-top 
source are given in Table|2] The LTS size was about 32 - 35 arc- 
sec in diameter and it was not changing significantly with time. 
The LTS had an elliptical shape with longer axis running along 
the tops of the loops. The loops heights and LTS altitude slowly 
increased. The altitude grew from 47 x 10^ to 60x 10^ km with an 
average velocity of 1 km s"'. The LTS spectrum could be fitted 
with one thermal component for all time intervals except for the 
three earliest ones for which we had to use a second component. 
The best fit was obtained for a thermal component with a power- 
law component (see Fig.|5] see also Conclusions for details). 

Estimated values of LTS heating rate are shown in Fig.|6] The 
actual value of the heating rate was about 0.1 -0.01 ergcm^-'s"'. 
Non-zero value of Eh is confirmed by diagnostic diagram (see 
Fig. Illl. The decay path of the flare runs along the QSS path 
i.e. there was continuous release of the energy and the rate of 
the release decreased very slowly which resulted in very slow 
cooling of LTS plasma. The characteristic time of temperature 
decrease was about 8 hours. 



4.2. The 22 August 2005 flare 

RHESSI light curves of the LDE2 are shown in Fig. |8] The flare 
emission could be analysed for ten intervals (see Table [1] and 
O. In the lower energy ranges 6-12 keV, RHESSI recorded the 
emission for 7 hours during the decay phase until two-hour gap 
in the satellite data. In the higher range 12-25 keV the emission 
was detectable for 3 hours after flare maximum. 
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Fig. 5. RHESSI X-ray spectrum for the loop-top source of the 30 
July 2005 flare recorded at 11:30 UT (horizontal bars, the bars 
widths corresponds to the energy bands widths). This spectrum 
was fitted using the thermal component, power-law non-thermal 
component and the spectral line complex at 6.7 keV (gray dashed 
lines). The total fitted model is marked by black line. 



Images of the flare and its evolution are shown in Fig. |9] 
(FIT 195 A images with LTSs observed by RHESSI overlaid on). 
The LDE2 was visible as an arcade of loops in EUV while in 
X-rays there was three loop-top source located just above the ar- 
cade. The first two sources (LTSl, LTS2) were visible for about 
4.5 hours until 6 UT. The third LTS appeared just before 6 UT 
and was present at least until the end of observations, i.e. for 
over two hours. The LTSl and LTS2 were not separated well 
enough to be analysed individually, therefore we treated them as 
one source denoted as LTS 1-2. 

Basic geometrical and physical properties of the loop-top 
sources are given in Table |3] The LTS 1-2 size grew from about 
30 arcsec to 50 arcsec in diameter while LTS3 size remained 
almost unchanged (40 arcsec, excluding the last observation 
point). The three loop-top sources were located along the tops 
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Fig. 4. SoHOfEVl 195 A images illustrating the decay phase of the 30 July 2005 flare. Contours show the emission in the 7-8 keV 
range observed with RHESSI. 
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Fig. 6. The 30 July 2005 flare - estimated range of the heating 
rate of the loop-top source. The real value of the heating must 
be contained between upper (upper symbols) and lower (lower 
symbols) limit of heating rate. 
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Fig. 7. GOES diagnostic diagram for the 30 July 2005 flare. The 
OFF and QSS paths of evolution are marked. The flare decay 
path runs along the QSS path. This confirms non-zero, slowly 
decreasing heating during decay phase of the flare. 
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Fig. 8. RHESSI light curves for the 22 August 2005 flare. 
Vertical lines mark the boundaries of the satellite night (N) and 
the SAA (S) periods (solid and dashed lines mark the beginning 
and the end of the periods, respectively). The light curves were 
shifted vertically for for clarity. The intensity is given in relative 
units. 



of the loops. The altitude of the LTSl-2 grew from 31 X 10^ to 
65 X 10^ km with an average velocity of 2 kms ' . The LTS3 was 
located higher but is moved upwards with velocity two times 
lower. The emission of each LTS could be well described by a 
single thermal component for the whole decay phase. 

The heating rate of LTS 1-2 and LTS3 are shown in Fig. [10] 
Eh was non-zero and slowly decreased with time as can be 
seen. Such temporal behaviour of the heating rate is confirmed 
by diagnostic diagram (see Fig. fTTT i. The decay branch of the 
flare runs along the QSS path or even flatter (due to altitude- 
density factor). The characteristic time of temperature decrease 
was about 6 and 7 hours for LTS 1-2 and LTS 3, respectively. 



4.3. The 25 January 2007 flare 

RHESSI light curves of the LDE3 are shown in Fig. [12] The flare 
emission could be analysed for twelve intervals (see Table [T] and 
|4]i. In the energy range 6- 12 keV, RHESSI recorded the emission 
for 6 hours during the decay phase until two-hour gap in the 
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Fig. 9. SoHO/EYT 195 A images illustrating the decay phase of the 22 August 2005 flare. Contours show the emission in the 7-8 keV 
range observed with RHESSI. 



Table 3. Basic parameters of the loop-top sources of the 22 
August 2005 flare. 
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Fig. 10. The 22 August 2005 flare - estimated range of the heat- 
ing rate of the LTSl-2 source (squares) and LTS3 source (di- 
amonds). The real value of the heating must be contained be- 
tween upper (upper symbols) and lower (lower symbols) limit 
of heating rate. 
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Fig. 11. GOES diagnostic diagram for the 22 August 2005 flare. 
The OFF and QSS paths of evolution are marked. The flare decay 
path runs along the QSS path (or even flatter). This confirms non- 
zero, slowly decreasing heating during decay phase of the flare. 



satellite data. In the higher range 12-25 keV the emission was 
detectable for 1 hour after the LDE3 maximum. 

Fig. [13] shows images of the LDE3 during its decay phase 
(BIT 195 A images with intensity isolines from RHESSI). The 
LDE3 was a limb arcade flare with only footpoints obscured by 
the solar disk (see FIT images). X-rays emission came from a 
source located just above FUV loops. In some energies and time 
intervals there were two LTSs, but there was too few such obser- 
vational points to analyse both LTSs separately and we decided 
to threat them as one source. 

Basic geometrical and physical properties of the loop-top 
source are given in Table |4] The LTS diameter changed from 
about 30 arcsec to 40 arcsec while its altitude grew from about 
30x 10^^ to 75 X 10^^ km with an average velocity of 2 kms ' . The 
emission of the LTS could be well described by a single thermal 
component for all time intervals. 

The heating rate of the LTS of LDF3 are shown in Fig. [14] 
Eh was non-zero and slowly decreased with time. This temporal 
behaviour of the heating rate is confirmed by diagnostic diagram 
(see Fig. [TsT l. The decay branch of the flare runs along the QSS 
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Fig. 12. RHESSI light curves for the 25 January 2007 flare. 
Vertical lines mark the boundaries of the satellite night (N) and 
the SAA (S) periods (solid and dashed lines mark the beginning 
and the end of the periods, respectively). The light curves were 
shifted vertically for for clarity. The intensity is given in relative 
units. 

Table 4. Basic parameters of the loop-top source of the 25 
January 2007 flare. 
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path or even flatter (due to altitude-density factor). The charac- 
teristic time of temperature decrease was very long i.e. 21 hours. 



5. Conclusions 

Previous results of the investigation of long-lasting HXR sources 
were based on the data giving rough information of physical pa- 
rameters of a source. The Yohkoh HXT L channel gave an in- 
tegrated flux in the wide energy band, 14 - 23 keV. Using the 
RHESSI data we are now able to analyse the HXR flux and 
images with energy resolution as good as 1 keV. We analysed 
RHESSI images and spectral data for the three selected long- 
duration flares which occurred on 30 July 2005, 22 August 2005 
and 25 January 2007. 

The flux for energy above 12 keV was detected for 6, 3 and 
1 hour during the decay phase of LDEl, LDE2 and LDE3, re- 
spectively. Below 12 keV emission was observed for at least 10, 
7 and 6 hours after maximum of LDEl, LDE2 and LDE3, re- 
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Fig. 14. The 25 January 2007 flare - estimated range of the heat- 
ing rate of the loop-top source. The real value of the heating must 
be contained between upper (upper symbols) and lower (lower 
symbols) limit of heating rate. 
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Fig. 15. GOES diagnostic diagram for the 25 January 2007 flare. 
The OFF and QSS paths of evolution are marked. The flare decay 
path runs along the QSS path (or even flatter). This confirms non- 
zero, slowly decreasing heating during decay phase of the flare. 



spectively. Such long-lasting high-energy emission is the first 
perceptible circumstantial evidence for the presence of heating 
processes during the decay phase. 

For all analysed LDEs and all time intervals we observed 
large hard X-ray loop-top sources. Each LTS observed was 
smooth, without any internal structure despite high angular res- 
olution of RHESSI (better than 10 arc sec). No LTS was visible 
in images reconstructed for a single narrower grid number 1-4 
(angular resolution 2.3 - 12 arcsec), so a large LTS cannot be 
explained as a superposition of smaller unresolved sub-sources 
unless separation the sub-sources is smaller than angular resolu- 
tion of the finest RHESSI grid (2.3 arcsec). 

The emission of the LTSs can be described as a single ther- 
mal component for almost all time intervals. Temperature of the 
component was quite high (more than 10 MK in the early part 
of the decay phase) and decreased very slowly. A characteristic 
time of the decrease was at least 6 hours which is several times 
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Fig. 13. SoHOfElT 195 A images illustrating the decay phase of the 25 January 2007 flare. Contours show the emission in the 
7-8 keV range observed with RHESSI. 



longer than a characteristic time of radiative cooling. A charac- 
teristic time of Spitzer conductive cooling is up to several min- 
utes. This is indirect but strong evidence that efficient heating of 
LTSs is needed to enable their long-lasting existence. 

This conclusion is confirmed by an energy balance calcu- 
lation. The actual value of Eh is larger than zero in all anal- 
ysed time intervals, implying that the energy release occurred 
many hours after the flares maximum. Total amount of ther- 
mal energy released during the decay phase is big and is about 
10^' erg cm"^ s"'. This is not slight amount of energy and it can 
be comparable with or even greater than the amount of energy 
released during the impulsive phase of a flare (see Jiang et al., 
l2006h . Non-zero, slowly decreasing heating during the decay 
phase was confirmed independently by the analysis of diagnostic 
diagrams. The decay evolution path runs along the QSS line or 
sometimes even flatter (due to altitude-density factor) for each 
of the three LDEs. 

Similar results concerning the energy release during the 
decay phase of LDEs were obtained in analysis based on 
Yohkoh/SXT data (e.g. Isobe et al., 2002; Kolomanski, 2007b). 
However, RHESSI data enable us to estimate much more reli- 
able value of plasma temperature - the SXT telescope had lim- 
ited sensitivity to hot plasma. Hence values of temperature and 
therefore heating rate obtained from RHESSI data are higher 
than those obtained earlier from the SXT data. A proper solar 
flare model should be able to explain these more severe results. 

In the case of LDEl for the three earliest time intervals the 
emission of LTS cannot be described by a single spectral compo- 
nent regardless of the simple structure of the source. The analysis 
of spatially resolved spectra of the LDEl showed that LTS emis- 
sion at that time consisted of two qualitatively different spectral 
components. The best fit to the observed spectra was obtained 
for a thermal component plus a power-law component with in- 
dex y - 9.7 and break energy equal 12 keV. A slightly worse fit 
is obtained for two thermal components. This result puts some 
requirements on solar flare models regardless of the nature of the 
second component. If both components are thermal, some mech- 
anism limiting fast mixing of hotter and cooler plasma in an LTS 
should exist. Moreover, the hotter component should be contin- 
uously heated by some mechanism. If one component is ther- 
mal and second is non-thermal, an acceleration process occur- 
ring within the hot (10 MK) and dense (10^ - 10"^ cm"^) region 
is needed. Both cases could be well explained by MHD turbu- 
lence in an LTS region. The turbulence produces large amount 
of small-scale reconnection regions where magnetic energy is 



continuously converted into thermal and non-thermal energy of 
plasma. Turbulence can also impede an escape of thermal en- 
ergy generated by these small-scale reconnection regions outside 
since in coronal conditions turbulent conductivity may be sev- 
eral orders of magnitude lower than classical thermal conduc- 
tivity parallel to magnetic field lines. The MHD turbule nce in- 
side LTSs wer e proposed b y several authors (e .g. Jakimiec et al.L 
119981: IJiang et al.L 120061: iKontar et all l201lh . A diff'erent ex- 
planation of two distinct thermal compon ents observed withi n 
a flare coronal source was proposed by ICaspi & LinI (|2010|) . 
Authors performed an analysis of short-duration flare and con- 
cluded that hotter (> 20 MK) component came directly from 
coronal reconnection region while cooler component ( < 20 MK) 
originates mainly from chromospheric evaporation. 

Here we reported three well observed events. In a future in- 
vestigation we will analyse a larger number of flares. Such analy- 
sis will give us statistical information on efficiency of the energy 
release in LDEs. 
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